The radiation— matter transition in Jordan— Brans— Dicke theory 



Andrew R. Liddle, Anupam Mazumdar and John D. Barrow 
Astronomy Centre, University of Sussex, Falmer, Brighton BN1 9QJ, U. K. 

(February 1, 2008) 

We study the transition from radiation domination to matter domination in Jordan-Brans-Dicke 
theory, in particular examining how the Hubble length at equality depends on the coupling parameter 
ui. We consider the prospects for using high-accuracy microwave anisotropy and large-scale structure 
data to constrain uj more strongly than by conventional solar system gravity experiments. 
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I. INTRODUCTION 

One of the most important epochs in the history of 
the Universe is the transition from radiation domination 
to matter domination. This transition alters the growth 
rate of density perturbations: during the radiation era 
perturbations well inside the horizon are nearly frozen 
but once matter domination commences, perturbations 
on all length scales are able to grow by gravitational in- 
stability. Consequently, the horizon scale at the time 
of matter-radiation equality is imprinted upon the spec- 
trum of density perturbations; indeed, in a flat cold dark 
matter cosmology it is the only length scale to appear in 
the perturbation spectrum. The nature of this transition 
has been well studied in the standard general relativistic 
cosmology, and plays a crucial role in calculations of the 
perturbation spectrum and associated microwave back- 
ground anisotropies by codes such as CMBFAST Q| . 

The success of general relativity as a description of 
our Universe allows us to evaluate the performance of 
rival theories of gravitation. So far, no weak-field or cos- 
mological observations disagree with the predictions of 
general relativity. The most important class of deviant 
theories are scalar-tensor gravity theories, of which the 
Jordan-Brans-Dicke (JBD) theory [|]|| is the simplest 
and best-studied generalization of general relativity. This 
theory leads to variations in the Newtonian gravitation 
"constant" G, and introduces a new coupling constant 
uj, with general relativity recovered in the limit Y/uj — > 0. 
The most robust constraint on uj, that it must exceed 
500, has been derived from timing experiments using the 
Viking space probe Q , and has stood for nearly 20 years 
now. Other constraints, such as those from nucleosynthe- 
sis ||), are comparable but more model dependent; the 
most detailed analysis Q gives only u> > 50. 

Within the next five to ten years, the advent of new 
microwave anisotropy satellites, MAP and Planck, and 
large galaxy-redshift surveys, 2df and the Sloan Digital 
Sky Survey (SDSS), promises to revolutionize our under- 
standing of cosmology by permitting the accurate deter- 
mination of a large number of cosmological parameters 
. So far, estimates of the accuracy of parameter 
estimation have only been made for cosmological param- 



eters, such as the Hubble constant Ho and the matter 
density f2 m , and for parameters describing the initial per- 
turbations, such as the spectral index n. However, such 
techniques can also in principle be extended to constrain 
parameters denning the underlying gravity theory, such 
as uj, which also influence the gravitational instability 
process. A full analysis of the viability of obtaining such 
constraints is an imposing task; all the perturbation for- 
malism employed to compute the present-day matter and 
radiation power spectra must be generalized to the the- 
ory of gravity under consideration, and the results then 
processed through the Fisher information matrix technol- 
ogy of Refs. In this paper, we assess whether such 
constraints might be competitive with existing bounds. 
We do this by studying the properties of the radiation- 
matter transition in JBD theory. We find that new cos- 
mological data sets may well give limits competitive with 
those obtained from weak-field solar system tests of gen- 
eral relativity. 



II. THE EQUATIONS 

The equations for a zero-curvature Friedmann universe 
are 
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where the Brans-Dicke coupling, lj, is a constant, a(t) is 
the cosmological scale factor, and is the Brans-Dicke 
field whose present value gives the observed gravitational 
coupling. Here, p and p are the energy density and pres- 
sure of the cosmic fluid, which has both matter (p m = 0) 
and radiation (p T — p r /3) components, so p = p m +p r and 
p = p T . Assuming negligible energy transfer between the 
fluids, they still obey the general- relativistic fluid conser- 
vation equations with solutions 
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where subscript '0' indicates the present value. We set 
a = 1. 

We shall assume a spatially-flat universe. It is impor- 
tant to realize that since the Brans-Dicke field appears in 
the Friedmann equation, the corresponding matter den- 
sity will differ from the general relativity 'critical' value, 
the correction being of order l/u), and this must be taken 
into account in our calculation. In fact, the present den- 
sity of matter corresponding to a flat universe is higher in 
JBD theory, as the Brans-Dicke field contributes a nega- 
tive effective energy density in the Friedmann equation. 

The present radiation density has two components. 
The photon energy density is dominated by the cos- 
mic microwave background, whose measured tempera- 



ture T-v = 2.728 ± 0.004 
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p lfi = 4.66 x 10~ 34 gem 
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However, the standard cosmology also contains relativis- 
tic neutrinos. These cannot be directly detected, and so 
their contribution must be fixed theoretically. Because 
neutrinos decouple before electron-positron annihilation, 
they are at a lower temperature than the photons by a 
factor {fljll (see e.g. Ref. fjj)), and, assuming three 
families of massless neutrinos, their total contribution is 
0.68 times that of the photons giving a present total den- 
sity in relativistic particles of 
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Accurate modelling of neutrino decoupling may affect 
this at the percent level, in a calculable way, as might 
contributions from unknown particles or gravitons decou- 
pling at much higher energies. The redshift of matter- 
radiation equality is given by 
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Our assumption of spatial flatness fixes the present mat- 
ter density, where one must be careful to take the contri- 
bution of the Brans-Dicke field into account when defin- 
ing the density parameter, fl. 



A. Analytic approximation 

When one of the fluid components dominates, the at- 
tractor solutions are well known jl3) . For radiation dom- 
ination, it is exactly the general relativity solution 
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while for matter domination there is a slow variation of 
the gravitational coupling described by the exact solu- 
tionis 



*In general, the JBD cosmologies have exact solutions which 
show that they are dominated by the $ field at early times and 
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Here, Eq. 

Newton gravitation constant to the $ field's value today, 
the relation being obtained by post-Newtonian expan- 
sion 101. Eq. (11) shows how the value of the matter 



density giving a spatially-flat universe is changed from 
the standard general relativistic expression by a finite to 
value. 

We are primarily interested in the Hubble radius at 
matter-radiation equality, which sets the characteris- 
tic scale of the bend in the fluctuation spectrum. It 
can be estimated analytically by assuming that the 
matter-dominated solution holds all the way to equal- 
ity; at equality, both the radiation and matter contribute 
equally to p, of course. We solve Eq. (|l|) with the help 
of Eqs. (H)-([l2|) to estimate the Hubble radius. We take 
only the first-order dependence on l/u> in our calculation, 
which yields 
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Using Eqs. (||) and ([li"|), we note that to O(l/co) we 
have 



(14) 



1 + z cq = 24000ft 2 ( 1 + — ] , 



where ft as usual is the present Hubble constant in units 
of 100 km s -1 Mpc , leading to 
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The leading term in Eq. ([L5|) is the general-relativistic 
limit, and is an exact result. Other terms in the equa- 
tion are the corrections accounting for the variation in $ 
between matter-radiation equality and the present. 



by the perfect fluid matter fields at late times. The general 
solutions approach the exact JBD vacuum solutions as t — > 
and the particular fluid solutions (given by Eqs. (p|)-(|l^) in 
the dust case) as t — > oo for a fiat universe. These attractor 
solutions also arise as simple exact solutions of Newtonian 
gravity with G <x t' n and a(t) oc t (2 ~ n)/3 , see Ref. @. 
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FIG. 1. A comparison of numerical and analytical solutions 
of the Hubble radius at equality for h — 0.5 and h — 0.7. The 
symbols correspond to the numerical results, and the smooth 
curves to the analytical estimates of Eq. ( |l5| ) valid in the large 
uj limit. 



B. Numerical approach 

The above result will be accurate at the large uj values 
that are of prime interest. However, to obtain the exact 
behaviour, we also tackle the problem numerically. 

In addition to the present values of a, a and we 
also require $ in order to specify a full solution of the 
JBD equations. In practice, this additional freedom is 
eliminated by the attractor behaviour during the long 
radiation-dominated phase. The attractor solution can 
be picked out in the initial conditions by requiring that 
$a approaches zero for a tending to zero at the initial 
singularity ||. Integrating Eq. once (taking advan- 
tage of the radiation contribution cancelling out), we ob- 
tain 
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and choosing the growing mode amounts to setting the 
constant C to be zero. 

Our numerical code establishes the smooth transition 
between the radiation- and matter-dominated solutions 
[Eqs. (0) and (|)]. We study the Hubble length at the 
transition for all positive values of uj and plot our numer- 
ical output in Fig. 1, for two different values of h. We 
compare our numerical result with the analytical esti- 
mate and confirm our analytical result at the asymptotic 
limit. At lower values of uj, the analytic approximation 
becomes inaccurate (due to the neglect of terms of or- 
der uj~ 2 and higher), and underestimates the change in 



III. PARAMETER ESTIMATION 

The shift in matter-radiation equality will lead to a 
shift of the maximum of the power spectrum. In the nor- 



mal general relativity situation, and assuming all dark 
matter is cold, this shift is governed by the combination 
of f£ m /i, where fi m is the matter density parameter, and 
this is often denoted T. Changing T gives a horizontal 
shift to the power spectrum P(k). Other parameters, 
such as the baryon fraction of the total matter density, 
also influence the location of the maximum, and some- 
times r is still used as an approximate description of this 
p"5[ , although future observations are expected to be so 
good that the baryons have to be included more accu- 
rately Jig ]. Other modifications to the cosmology, such 
as inclusion of a hot dark matter component, cannot be 
modelled solely by a shift in the power spectrum. 

At the current 95% lower limit of uj = 500, the shift 
in matter-radiation equality will be at the one percent 
level. The parameter uj can also be probed via its affect 
on the growth rate of perturbations although, because 
of biasing, galaxy surveys have problems in constraining 
the overall normalization. Note that the lack of a char- 
acteristic scale during matter domination means that the 
altered growth rate will not change the shape of the spec- 
trum; this result holds in JBD theory as well as in general 
relativity. 

Measurement of the shape of the power spectrum 
from galaxy surveys has been considered in Refs. ^,|). 
Tegmark |8|] introduces a phenomenological parameter r\ 
which shifts the power spectrum horizontally, exactly the 
effect we are interested in, and considers how accurately 
it can be measured by the SDSS. The accuracy depends 
on the shortest scales considered; one cannot go too far 
without worrying about non-linear clustering and biasing 
effects. If we take Tegmark's estimates evaluated at the 
scale currently going non-linear (/c non -iin — 0.1 /iMpc" 1 ), 
we can read off the anticipated errors under two assump- 
tions. If all the model parameters are to be determined 
from the SDSS alone, A77/77 ~ 0.1. If, on the other hand, 
it is assumed that all parameters are already fixed except 
77, then A77/77 ~ 0.02. This is entering at about the re- 
quired level. At least in principle, data on scales in the 
non-linear regime can improve this further if their theo- 
retical interpretation is deemed sufficiently unambiguous. 

The case where other parameters can be considered as 
fixed might well apply once one takes microwave back- 
ground anisotropics into account, and uses them to com- 
pute the other cosmological parameters. Indeed, if gen- 
eral relativity is correct, then the Planck satellite alone 
will already be able to measure 77 at the percent level 
and then galaxy surveys can be used to improve the 
estimate further. However, there are many potential de- 
generacies: the horizon size at equality is changed by 
altering Q m , or h, or increasing the number of massless 
species (perhaps even by including a possible thermal 
graviton background) as well as by introducing finite uj. 
On the other hand, the degeneracy may be broken by 
the different growth rate of perturbations in JBD the- 
ory. The evolution of dust density perturbations (where 
5 = 5p/ p <C 1) for all wavelengths is determined by the 
solution of IEtLCoI 
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Hence, the growing mode solution of Eq. ( [l7|) for the 
background Universe of Eqs. (0) and <M) is given by 
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This gives a normalization shift at the few percent level 
for u> ~ 500, most of the effect being after the microwave 
anisotropies have been generated, since the redshift of 
decoupling is close to that of equality. 

It appears that there is a reasonable prospect that up- 
coming precision observations can impose a limit on (or 
make a detection of) finite u> values at a level competi- 
tive with the post-Newtonian bounds, although the use 
of galaxy surveys for thispurpose is subject to a number 
of caveats given in Ref . || . While we have presented the 
case for arriving at new limits on the constant value of u> 
characterising JBD theories, the basic technique can be 
extended to constrain the value of a non-constant uj($) 
defining cosmological solutions to a more general scalar- 
tensor gravity theory, as discussed in Ref. Jl8f . 



IV. CONCLUSIONS 

We have studied the matter-radiation transition in the 
JBD theory, both numerically and analytically. The shift 
in the epoch of matter-radiation equality will influence 
the shape of the density perturbation spectrum, and it 
appears that precision microwave anisotropy measure- 
ments and large galaxy-redshift surveys may in the future 
be able to impose limits on to competitive with existing 
solar system bounds. 

However, it seems unlikely that a very substantial im- 
provement will be possible. It may therefore be best to 
wait to see whether the high quality of data promised 
is actually delivered before embarking on the substantial 
undertaking of generalizing general-relativistic results to 
carry out a proper estimate of the likely observational 
limits. If all goes well, use of the data to constrain pa- 
rameters of the gravitational theory will be a worthwhile 
endeavour and an unexpected bonus from future high- 
precision observational studies of galaxies and the mi- 
crowave background. 
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